Many body effects in finite metallic carbon nanotubes. 
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The non homogeneity of the charge distribution in a carbon nanotube leads to the formation of an 
excitonic resonance, in a similar way to the one observed in X-ray absorption in metals. As a result, 
a positive anomaly at low bias appears in the tunnelling density of states. This effect depends on 
the screening of the electron-electron interactions by metallic gates, and it modifies the coupling of 
the nanotube to normal and superconducting electrodes. 

PACS numbers: 73.63.Fg , 73.23.Hk 
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Introduction. The addition of single electrons to fi- 
nite carbon nanotubes induces measurable effects, re- 
lated to the finite spacing between the energy levels 
and to the electrostatic energyassociated to the electron 
charged H II llJUzl B 13 The standard model for 
Coulomb blockade[l(J |Tl| assumes that the electrostatic 
potential inside the system is raised by the addition of an 
individual electron, preventing tunnelling unless different 
charge states are degenerate. The potential is supposed 
to be constant throughout the sample. When it is not the 
case, non equilibrium effects may occur [T^ll3| . which are 
related to the Fermi edge singularities associated to the 
sudden ejection of a core electron in a metal 14]. One 
dimensional systems, like the nanotubes, are good can- 
didates for observing these effects, as screening is sup- 
pressed, and the electrostatic potential inside them can 
be modulated by metallic gates in their proximity. Other 
effects related to modulations in the electrostatic poten- 
tial were considered infill. 

We calculate in this paper the changes in the tunnelling 
density of states due to the non homogeneity of the elec- 
trostatic potential, for different possible experimental se- 
tups. The effects discussed here look similar, but are dif- 
ferent from the Luttinger liquid features expected near a 
contact fl6l [ill [ill. fUS f . which have also been observed in 
nanotubes |20l |21| . 

The next section describes the model to be studied. 
The main results are presented next. Then, we analyze 
how to incorporate Luttinger liquid effects. The expected 
behavior at energies comparable to the spacing between 
individual electronic levels is discussed next. The paper 
concludes with a section on possible experimental conse- 
quences of the effects analyzed here. 

The model. We consider only a nearest neighbor hop- 
ping, t, between n orbitals at the C atoms. Each sub- 
band in a zigzag nanotube can be modelled by a one 
dimensional tight binding hamiltonian with two sites per 
unit cell and two different hoppings, t and t(l + e 1 *") — 
2te i< f>"/ 2 cos(0„/2), where qb n = (2wn)/N, n = 1, • • • , N- 
1 is the transverse momentum associated to the subband, 
and N is the number of C atoms in a transverse section 
of the nanotube. 

In a metallic nanotube, N = 3M, where n is an in- 
teger. There are two gapless bands, characterized by 



FIG. 1: Sketch of one of the geometries considered in order 
to calculate the charge distribution in the nanotube. 



n = ±M. Other bands with n ^ M have a gap, 
A„ = f |1±2 cos(0„/2)|. Thus, the electronic states of the 
two bands which cross the Fermi level can be described 
by the simple hamiltonian: 
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The Fermi velocity can be written as vp — 2ta, where a is 
the lattice constant. We consider a finite nanotube with 
L unit cells. When an electron hops into the nanotube its 
charge will be distributed throughout its length. In a first 
approximation, this charge can be calculated using the 
Hartree or Hartree-Fock approximation. As the electro- 
static potential is weakly screened in a one dimensional 
geometry, the charge distribution can be influenced by 
metallic gates in the vicinity. We first consider the setup 
depicted in Fig.pP, where the electron is emitted from 
a metallic electrode at one end of the nanotube. More 
complicated geometries will be studied later. 

We consider the two spin degenerate bands which cross 
the Fermi level. The polarization of the remaining bands 
by the potential associated to the charge of a single elec- 
tron will be small, and it can be treated perturbatively. 
The radius of the nanotube acts as a short distance cut- 
off of the Coulomb potential. We describe the effects 
of the metallic electrode by an image charge induced by 
the physical charge on the nanotube. We approximate 
the electrostatic interaction between electrons located at 
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FIG. 2: Electrostatic potential and charge (inset) at the edge 
of a nanotube with an additional electron and L = 1024 unit 
cells. The parameters used are e 2 /«F = 5.4, R/a — 3 and 
d/a = l. 



unit cells n\ and ni as measured from the gate, as: 



Vtot(«i,n 2 ) = V c i(ni - n 2 ) - V e \ ni + n 2 + - 

a 



V"d -ni - ri2 + - + 14i(-ni + n 2 ) 



V d (n) 



Vn 2 a 2 + R 2 



(2) 



where d is equal to twice the distance between the nan- 
otube and the electrode, and R is the radius of the nan- 
otube. 

Results. The electrostatic potential and the induced 
charge when the number of electrons in one of the four 
subbands at the Fermi level is one above half filling 
(which is taken as the neutral situation) is shown in 
Fig.0. 

The calculations have been done using the Hartree 
approximation. The results are significantly changed 
if exchange is included. We find enhanced Friedel 
oscillations 22], and a sizable gap pinned at the Fermi 
level. Some of these effects are due to changes in bulk 
properties which are unrelated to the addition of elec- 
trons. In the following, we present results obtained 
within the Hartree approximation, where the features as- 
sociated to single charges are easier to isolate. 

A non negligible fraction of the charge is localized by 
its image near the electrode. This leads to a reduction of 
the repulsive electrostatic potential in that region. This 
effect is relatively small compared with the gaps of the 
higher lying subbands of a small carbon nanotube, so 
that the assumption of neglecting their polarization is 
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justified. We also find substantial Friedel oscillations, 
which have periodicity two, as fcp = n/a. 

The charge distribution depends strongly on the loca- 
tion of the external electrodes. Fig. [3] compares the po- 
tential calculated previously with the one obtained when 
there are two symmetrically placed electrodes at each 
end of the nanotube, and in the absence of electrodes, all 
other parameters being the same. 
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FIG. 3: Electrostatic potential induced by a single elec- 
tron when there is one electrode (full line), two symmetri- 
cally placed electrodes (dashed line), and no electrodes (dash- 
dotted line). The inset shows the non equilibrium density of 
states for the same three cases. 



The effects associated with the non homogeneous 
charge distribution after the injection of a single elec- 
tron are described, within the Hartree approximation 
used above, as a change of all electronic levels. Taking 
this effect into account, the electronic density of states 
is changed as: i) The overlap between the initial and fi- 
nal eigenstates is reduced. In the limit of vanishing level 
spacing, this leads to the orthogonality catastrophe in- 
duced by the sudden switching of a local potential [23|. 
ii) The potential shifts the electrons towards the elec- 
trode, enhancing the density of states at the Fermi level 
in its vicinity. This is the so called excitonic effect [bl Izif 
and it opposes the orthogonality catastrophe. 

The final local density of states, in the limit of vanish- 
ing level splitting, goes as D(uo) oc cj^ 1- " 5 ) 2-1 ]^!, where 
5 > 0. Hence, for sufficiently small perturbations, 6 < 1, 
the tunnelling density of states is enhanced at the ends of 
the nanotube. The effective density of states can be ob- 
tained by assuming that the hamiltonian after the charg- 
ing of the nanotube is modified from an initial expression 
Ho to Hf = Ho + SV SC , where 5V SC is the modification of 
the Hartree-Fock potential induced by the extra electron. 
Assuming that the tunnelling takes place at position 1 
and subband 1 near the edge of the nanotube, the tun- 
nelling density of states at zero temperature is: 



GH = II I Ai<°> N \°> N )o,jf x E I f^ n > N + ^olO, N) ,i 8{uo - e n + e ) 



(3) 



3 



1 ' 
□ 




0.05 0.10 0.15 

Energy ( v p / a ) 



FIG. 4: Effective tunneling density of states at the end of a 
nanowire with L = 1024 unit cells, using the same parameters 
as in Fig. [2| . The dotted line is the density of states in the 
absence of final state effects. The inset gives the density of 
states over the entire bandwidth of the conduction band. 



where the states \n, N)^j),j are eigenstates of H f with 
N electrons, and j = 1, • • • 4 is the band index. As these 
states are Slater determinants, the calculation is reduced 
to a sum of determinants of overlaps between one particle 
wavefunctions. 

Results obtained for the injection in a neutral nan- 
otube with L — 1024 unit cells and the parameters which 
lead to the electrostatic potentials shown in Fig.p] are 
shown in the inset of Fig.jSj, and, in more detail, in 
Fig- IH i where they are also compared with the results 
obtained neglecting <5V^, C - 

Luttinger liquid effects. The calculation discussed 
above takes into account, in an approximate way, the in- 
teraction corrections expected at the boundary of a Lut- 
tinger liquid. The anomalous exponent in the density 
of states can be estimated using first order perturbation 
theory, and it is finite within Hartree-Fock theory |2^|. 
Note that the calculations also describe approximately 
the crossover to a Fermi liquid behavior at sufficiently 
high energies plj. 

For sufficiently long wires, we can use the bosoniza- 
tion approach. The calculation of the final state effects 
leading to eq.© can be extended without too much dif- 
ficulty to the case of a Luttinger liquid with forward in- 
teractions. The initial and final bosonized hamiltonians 
arc 261: 



Ho = i / dx4irv F n 2 (x) + ^ld x $( 

2 Jx>0 47r 



Hi — Hq- 



8tt 2 



2tt 



dx \ dyd x $(x)Vtot(x,y)d y $>(y) 



(4) 



x>0 Jy>0 



dx6V BC (x)d x $(x) 



x>0 



The hamiltonians Ho and Hi differ by a term which is 
linear in the boson fields, and there is a simple unitary 
transformation, hi which transforms one into the other. 
The calcuation of the density of states can be written as 



the Fourier transform of the expression: 

G{t) = {Q\ijj{x = 0)e mt ^(x = 0)|0) 

= (0\ip(x = Q)U- l e iHot U^{x = 0)|0) (5) 

where |0) is the ground state of Ho, and we are assuming 
that the end of the nanotube is at position x = 0. 

The hamiltonians in eq.Q are quadratic in the boson 
coordinates, and the expression in eq.JSJ) can be calcu- 
lated by normal ordering the operators, which depend 
exponentially on the bosonic degrees of freedom. If wc 
neglect the breakdown of translational symmetry induced 
by the gates, and the long range effects of the Coulomb 
potential, the tunnelling density of states becomes: 



D(u>) oc \u\ g 



(6) 



where g is the parameter which describes 
the Luttinger liquid properties, and SV = 
(27r) _1 limfc^o f Q °° e lkx 5V sc (x)dx. For repulsive in- 
teractions, g < 1, and Luttinger liquid effects tend to 
suppress the positive anomaly studied here. 

In the presence of gates, the electron-electron poten- 
tial is not translationally invariant. The screening effects 
of the gates reduce locally the interactions, which are un- 
screened away from the gate. A related situation was con- 
sidered in|27j , a Luttinger liquid connected to non inter- 
acting gates. The quadratic hamiltonian, eq.Q, which is 
defined on a half line, x, y > 0, can be solved by defining 
a related problem on the entire axispsj. We can estimate 
in a simple way the value of the parameter g at the edge 
of the nanotube from the energy associated to a fluctua- 
tion of the charge density on a scale I < L in that region. 
This estimate, in the thermodynamic limit, leads to the 
charge compressibility and to the bulk value of g. In the 
absence of screening by gates, the inverse compressibility 
diverges logarithmically with the size of the system. It 
is easy to show that, in the geometry depicted in Fig. [I], 
the energy associated to a fluctuation of size I is given by 
a term which diverges logarithmically as I — > oo minus 
a term which tends to a constant. The second term de- 
scribes the screening effects of the gate. Thus, for large 
nanotubes and at sufficiently low energies, bulk effects 
determine the value of g, and the tunnelling density of 
states will be given by eq.®. At high energies, on the 
other hand, the screening of the electron-electron inter- 
actions by the gate makes g — > 1 and the non equilibrium 
effects described here will dominate. The crossover be- 
tween these two regimes depends on the limiting value of 
the interaction, which, in turn, is influenced by screening 
by other parts of the system (see below). 

Single electron properties. At sufficiently low ener- 
gies or temperatures, a finite nanotube behaves like a 
quantum dot [J H H H H HI 113-For the parame- 
ters used here, e 2 /vf — 5.4, the level splitting and the 
charging energy are comparable. The Hartree calcula- 
tion described above allows us to estimate the charging 
energy, Ec — En+i + En-i — 2 .En- The Coulomb po- 
tential considered here, in the absence of gates, leads 
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to Eq ~ e 2 L~ 1 \og(L/R). The screening by the vertical 
gates, as in Fig.pQ, reduces this value by a quantity which 
scales as L~ l . For the parameters used in Fig. [21 we find 
a level spacing (in units of v-p/d) of AE = 0.00306, and 
a charging energy Eq, = 0.0196. 

The final state effects considered here manifest them- 
selves, at the scale of the level spacing, as a dependence of 
the tunnelling amplitude on the charge state. This effect 
can be included by adding assisted hopping terms to the 
hamiltonian pni 132. l33| (see also (34) 1. This effect, in our 
case, is rather small. Using the parameters in Fig. [2] , the 
value of the square of the renormalized creation operator 
(see eq.Q) at zero energy is 0.002616 for the transition 
N -> N + 1 and 0.002612 for N + 1 -> N + 2 (note that 
these transitions involve the same electronic level). The 
smallness of the final state effects on single particle levels 
shows that the enhancement of the effective tunnelling of 
states is a collective phenomenon, which arises from the 
cumulative changes in many levels near the Fermi energy. 

Conclusions. We have analyzed, semi-quantitatively, 
the appearance of a low energy peak in the tunnelling 
of states in nanotubes due to the inhomogeneous poten- 
tial induced by the charge of a single electron. This zero 
bias anomaly is closely related to the excitonic resonance 
which appears in X-ray absorption in metals Qjjij , The 
formation of this resonance enhances the tunnelling den- 



sity of states, although it is induced by the Coulomb in- 
teraction. The strength of this anomaly depends on the 
screening of the nanotube by metallic contacts and gates, 
and, presumably, it can be tuned experimentally. 

The strength of the enhancement of the tunnelling den- 
sity of states can be suppressed by Luttinger liquid ef- 
fects, for repulsive electron electron interactions. Lut- 
tinger liquid effects are reduced when the interaction is 
screened, which, on the other hand, favors the existence 
of the excitonic resonance. The effects discussed here 
should also be present in multiwall nanotubes. In a sys- 
tem with many bands at the Fermi level, only one con- 
tributes to the resonance, while the others induce an or- 
thogonality catastrophe (note that we have considered 
four bands here, see eq.©). 

The presence of an excitonic resonance makes the 
transport properties of the nanotube similar to those of 
Luttinger liquids with attractive interactions. It can lead 
to the enhancement of the proximity effect in nanotubes 
attached to superconducting electrodes [3^ |3(| |37]], and 
increase, in general, the tendency towards superconduc- 
tivity in these systems [3^. 
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